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Abstract 

We calculate in heavy baryon chiral perturbation theory the local NN-potentials gener- 
ated by the two-pion exchange diagrams at two- loop order. We give explicit expressions for 
the mass-spectra (or imaginary parts) of the corresponding isoscalar and isovector central, 
spin-spin and tensor NN-amplitudes. We find from two-loop two-pion exchange a sizeable 
isoscalar central repulsion which amounts to 62.3 MeV at r = 1.0 fm. There is a similarly 
strong isovector central attraction which however originates mainly from the third order 
low energy constants dj entering the chiral 7riV-scattering amplitude. We also evaluate the 
one-loop 27r-exchange diagram with two second order chiral 7r7riViV-vertices proportional 
to the low energy constants c\$ 34 as well as the first relativistic 1/M-correction to the 
27r-exchange diagrams with one such vertex. The diagrammatic results presented here are 
relevant components of the chiral NN-potential at next-to-next-to-next-to-leading order. 
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t^- ' Over the past years effective field theory methods have been successfully applied to the two- 

nucleon system at low and intermediate energies [[!], 0, |3], |], [5|, ]E|]. The idea of constructing 
the NN-potential from effective field theory was put forward by Weinberg [|ij and this was first 
taken up by van Kolck and collaborators |2j who used "old-fashioned" time-ordered perturbation 
theory. Later, the systematic method of unitary transformations was employed by Epelbaum, 
Glockle and Meifiner Q to construct an energy-independent NN-potential from the effective 
chiral pion-nucleon Lagrangian at next-to-leading order. Based on one- and two-pion exchange 
and nine adjustable NN-contact interactions which contribute only to S- and P-waves a good 
description of the deuteron properties as well as the NN phase-shifts and mixing angles below 
Tiab = 300 MeV was found in that framework |3J . Also recently, the elastic proton-proton 
scattering data base below 350 MeV laboratory kinetic energy (consisting of 1951 data points) 
has been analyzed in terms of l7r-exchange and chiral 27r-exchange at next-to-next-to- leading 
order in ref.0. The resulting good x 2 /dof < 1 constitutes a convincing proof for the presence 
of the chiral 27r-exchange in the long-range proton-proton strong interaction. It was concluded 
in ref . @ that l7r-exchange together with chiral 27r-exchange gives a very good NN-force at least 
as far inwards as r = 1.4 fm internucleon distance. All shorter range components of the NN- 
interaction have been effectively parametrized in ref.Q by 23 boundary condition parameters. 

At present, there is much interest in extending the calculations of the two-nucleon system 
|| (and also the analysis of the NN-data base) to one higher order in the chiral expansion. 
This requires the full knowledge of the next-to-next-to-next-to-leading order (N 3 LO) chiral NN- 
potential. In momentum space this corresponds to terms in the NN T-matrix which are of fourth 
order in small external momenta and the pion mass, denoted generically by 0(Q A ). In particular, 
all two-loop diagrams of the process NN — > NN with leading order vertices contribute at 
that order, 0(Q A ). In three recent publications |J we have calculated completely the static 
NN-potentials generated by the (two-loop) 37r-exchange diagrams with all possible interaction 
vertices taken from the leading-order effective chiral 7riV-Lagrangian. In these works we have 
evaluated the mass-spectra or imaginary parts from which one can also easily reconstruct the 
momentum space NN-amplitudes in the form of a subtracted dispersion relation. The appearing 
subtraction constants can be absorbed in the strengths of some local NN-contact interactions 



which will be treated as adjustable parameters in a fit to NN-phase shifts etc. The purpose of 
this work is to present analogous results for the two-loop contributions to the chiral 27r-exchange 
which also belong to the chiral NN-potential at order 0(Q 4 ). Furthermore, we evaluate several 
new contributions at the same order 0(Q A ). These are the one-loop 27r-exchange diagram with 
two second order chiral 7T7riViV-vertices proportional to the low energy constants 01,2,3,4 and the 
first relativistic 1/M-corrections to the 27r-exchange diagrams with one such vertex proportional 

to Ci, 2,3,4- 

Let us first give some basic definitions in order to fix our notation. In the static limit and 
considering only irreducible diagrams the on-shell NN T-matrix takes the following form 
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where q = \q\ denotes the momentum transfer between the initial and final-state nucleon. The 
subscripts C, S and T refer to the central, spin-spin and tensor components, each of which 
occurs in an isoscalar (Vc,s,t) and an isovector version (Wc\s,t)- As indicated, the (real) NN- 
amplitudes Vc(q), ■ ■ ■ , Wr(q) depend only on the momentum transfer q in the static limit M — > 
00, where M denotes the nucleon mass. We are interested here only in non-polynomial or 
finite range terms. For this purpose it is sufficient to calculate the imaginary parts of the 
NN-amplitudes, ImVc,s,T{ifi) and lmWc,s,T{i/J'), which result from analytical continuation to 
time-like momentum transfer q = ifi — + with fi > 2m n . These imaginary parts are also the 
mass-spectra entering a representation of the local coordinate-space potentials in the form of a 
continuous superposition of Yukawa-functions, 
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The isoscalar central, spin-spin and tensor potentials, denoted here by Vc,s,t{t), are as usual 
those ones which are accompanied by the operators 1, &i -a 2 and 3 B\ -r a 2 -r — d\ -a 2 , respectively. 
For the isovector potentials Wc srij') a completely analogous representation holds, of course. 






Fig.l: In -exchange diagrams at two-loop order. The grey disc symbolizes all one-loop diagrams 
of elastic nN -scattering. The combinatoric factor of the first two diagrams is 1/2. Diagrams for 
which the role of both nucleons is interchanged lead to the same NN-potential. 



The two-loop 27r-exchange diagrams are symbolically represented in Fig. 1. The grey disc 
should be interpreted such that it includes all one-loop diagrams of elastic 7riV-scattering. The 
imaginary parts entering eqs. (2,3,4) can be calculated from these 27r-exchange diagrams as in- 
tegrals of the A^A^ — > 2tt — > A^A^ transition amplitudes over the Lorentz-invariant 27r-phase 



space making use of the Cutkosky cutting rule. In the present case these transition ampli- 
tudes are products of tree-level 7riV-amplitudes and one-loop 7riV-amplitudes. The pertinent 
27r-phase space integral is most conveniently performed in the nn center-of-mass frame where 
it becomes proportional to a simple angular integral f_±dx . In the center-of-mass frame the 
on-shell pion four-momenta read k\ 2 = (/i/2, ±k) and as explained in ref.|J the four-momenta 
of the external nucleons can be chosen as P" 2 — (i/^/2,pw ) with v a real unit vector and 

iM + 0(M~ 1 ) in the heavy nucleon limit. As an extensive check on these 



p = ^74 - M 2 

techniques we have recalculated the imaginary parts of the one-loop 27r-exchange diagrams (con- 
sidering graphs with no, single and double A-excitations, etc.) and we reproduced indeed exactly 
the results of refs. 0, |] ina rather short calculation. 

Let us now turn to the results for the 27r-exchange at two-loop order. From the last two 
diagrams in Fig. 1 we obtain the following imaginary part for the isoscalar central NN-amplitude, 
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Note that the expression in the curly bracket of eq.(5) is proportional to the isospin-even non- 
spin-flip one-loop 7riV-amplitude Reg^ op (0, fi 2 ) (see appendix A in ref.[|10|). Modulo the last 
term proportional to g\ this is precisely the result of the so-called correlated 27r-exchange given 
in section 4 of ref. ||. Of very similar structure is the imaginary part of the isovector spin-spin 
(and tensor) NN-amplitude, 
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Again, the expression in the curly bracket of eq.(6) is proportional to the isospin-odd spin-flip 
one-loop 7riV-amplitude Re/i z ~ op (0,/i 2 ) |I0| . Besides the last 5^4-term, eq.(6) agrees with result 
of correlated 27r-exchange given in ref.0. Next, we come to the imaginary parts of the isoscalar 
spin-spin and tensor NN-amplitudes. We obtain for them the following representation, 
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the pion center-of-mass momentum. Modulo the low-energy constant 



with k = 

^izi(A) — d[ 5 (\) introduced in ref. fill the expression in the curly bracket of eq.(7) is proportional 
to the isospin-even spin-flip one-loop 7riV-amplitude Re [h^ (ui, fj?)/uj] evaluated at uj — ikx. 
All four diagrams in Fig. 1 contribute to the isovector central NN-amplitude Wc{q) and we find 
for its imaginary part the following representation, 

2k n 



ImW c (ifi) 



M^m 



dx 



,2^2 



^(2m 2 -/i 2 ) + 2(^-l)A; 2 x 



2^2 



*jr\j X 



21n^ 
A 



1 



+Qkx\/m 2 + k 2 x 2 In ■ 



kx + \ ml + k 2 x 2 



m„ 



x 



+ 



5 m^ 
In — — 

6 A 



+ 



mt 



k 2 x 2 



(i + 



Amid 



2g 2 A ) 



// 2 (1 + 5g\ 



ml 
k 2 x 2 
k 



+ g A (fi 2 - 2k 2 x 2 - 2m\ 
3/2 kx + Jm 2 + k 2 x 2 ' 



ln 



m„ 



In 



// 



H + 2k 

2m 1T 



-^(l + 5^)ln^ 



^ 



+ ^(5 + 13^) -2m 2 (1 + 2^) 



+967T 2 / 2 (2m 2 - /i 2 )K(A) + 4(A)) - 2k 2 x l d\{\) + 4m 2 rf r 5 (A) 



Again, modulo the low energy constants d r A\) the expression in the curly bracket of eq.(8) is pro- 
portional to the isospin-odd non-spin-flip one-loop 7riV-amplitude Re [gf oop (ui, {i 2 )/u} evaluated 
at <jj = i kx. Furthermore, when restricting the second factor in eq.(8) to the terms in fourth and 
fifth line one recovers the result of correlated 27r-exchange given in section 4 of ref . . The last 
line in eq.(8) gives the contribution of the third order low energy constants d\ 2 3 5 (A) introduced 
in ref. |TT|. Such counterterms are necessary in order to absorb the divergences generated by the 
one-loop graphs of elastic 7riV-scattering. In fact, the chiral logarithms h^m^/A) and the scale 
dependent low energy constants d T j(\) in eqs.(7,8) combine to the (scale independent) barred 
quantities cL- defined in ref. |ll], or equivalently dj = dy 
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This completes the presentation 
of analytical results for the two-loop 27r-exchange NN-interaction. Modulo polynomials the mo- 
mentum space NN-amplitudes Vc(q), ■ ■ ■ , Wt{q) can be obtained in the form of a subtracted 
dispersion relation 
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For the isovector amplitudes Wc,s,t{^) a completely analogous dispersion relation representation 
holds, of course. 
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Tab.l: Numerical values of the local NN-potentials generated by two-loop chiral 2 r n-exchange 
versus the nucleon distance r. The units of these potentials are MeV. 

In Table 1, we present numerical values for the coordinate space potentials generated by chiral 
27r-exchange at two-loop order inserting eqs. (5,6,7,8) into eqs. (2,3,4). We use the parameters 
f n = 92.4 MeV, m^ = 138 MeV (average pion mass) and g^ = 1-3. Via the Goldberger-Treiman 
relation this corresponds to a strong 7riViV-coupling constant of g^N = 9aM/ 'f„ = 13.2 which is 
in agreement with present empirical determinations of g^N from 7riV-dispersion relation analyses 
fl2"f . For the third order low energy constants dj we use the average values of three fits to tcN- 
phase shift solutions given in ref.|ll[|: d\ + d 2 — 3.0 GeV -2 , d 3 = — 3.0GeV -2 , J 5 = 0.1 GeV~ 2 , 
du — G?i5 = —5. 7 GeV -2 . One observes from Table 1 that the isoscalar central potential Vc{r) 
generated by two-loop 27r-exchange is sizeable and repulsive. This repulsive potential is almost a 
factor 100 larger than the one from the so-called correlated 27r-exchange investigated in section 
4 of ref. . The origin of this strong enhancement is the additional and obviously dominant g\- 
term in the curly bracket of eq.(5) which stems from all those two- loop diagrams which were not 
considered in ref.0. Similar features hold for the attractive isovector spin-spin and the repulsive 
isovector tensor potentials Ws,t[t) which are however typically a factor 10 smaller than isoscalar 
central potential Voir). At first sight, the attractive isovector central potential Wc{r) from two- 
loop 27r-exchange seems to be even larger in magnitude. However, in this case the effect comes 
mainly from the low energy constants dj. For example, if one omits the d r AX) in eq.(8) and sets 
the scale A = m u = 782 MeV in the chiral logarithm ln(m„./A) the isoscalar central potential 
Wc{r) gets reduced by more than one order of magnitude. The same phenomenon is observed for 
the isoscalar spin-spin and tensor potentials Vs,T(r) which is also dominated by the contribution 
from the low-energy constant d\& — <i 15 . In the light of these findings one can expect that the 



effect of including the two-loop 27r-exchange in NN-phase shift calculations is not just a small 
correction. Of course, a firm conclusion about the size of such corrections can only be drawn 
from the complete N 3 LO chiral NN-potential. 
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Fig. 2: The heavy dot symbolizes the second order chiral irnNN -contact vertex proportional to 
the low energy constants 01,2,3,4- The 2n-exchange diagrams for which the role of both nucleons 
is interchanged are not shown. 

Furthermore, we like to present here results for another class of new contributions to the NN- 
potential at order 0(Q 4 ). This includes the one-loop 27r-exchange diagram with two second order 
chiral 7T7riViV-vertices proportional to the low energy constants 01,2,3,4 and the first relativistic 
1/M-corrections to the 27r-exchange diagrams with one such vertex proportional to 01,2,3,4. The 
corresponding graphs are shown in Fig. 2. The first "bubble" or "football" diagram can be 



straightforwardly evaluated in the heavy baryon formalism and one finds [13 
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In eqs. (10,11) and in the following we omit purely polynomial terms of the form const q 4 + 
const q 2 + const in the central NN-amplitudes and const q 2 + const in the tensor and spin-orbit 
NN-amplitudes. The third diagram in Fig. 2 with one ci,2,3,4-vertex starts to contribute at order 
0(Q 3 ) (see section 4 in ref. ||) and the next higher order contribution at 0(Q 4 ) arises from the 
first relativistic 1/M-corrections to the vertices and nucleon propagator. In order to be specific 



we employ here the minimal relativistic second order 7riV-Lagrangian C^j^ 1 defined in ref. |TT] 



A convenient way to obtain the 1/M-corrections is to use relativistic nucleon propagators and 
interaction vertices and to perform (in the NN center-of-mass frame) the 1/M-expansion inside 
the pion-loop integral. Consider now the second diagram in Fig. 2 (and its mirror partner). 
Since the left hand side Tomozawa- Weinberg vertex is of isovector nature only the isovectorial 
C4-vertex can contribute on the right hand side. One finds modulo polynomials, 
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As a genuine relativistic effect we encounter here a spin-orbit NN-amplitude Wso(q)- According 
to the convention of ref. 0, [| the spin-orbit amplitudes are accompanied in the NN T-matrix 
by the (momentum space) spin-orbit operator i{a\ + 0^2) • (qx p), where p denotes the nucleon 



center-of-mass momentum. The C4-vertex in the third diagram of Fig. 2 (and its mirror partner) 
generates further isovector NN-amplitudes, 

Wr(q) = -^W s {q) = ^^ (16^ + 7q 2 )L( q ) , (14) 

Wsob) = ~W c (q) = "^f^ « + tf)L(<l) ■ (15) 

Finally, the isoscalar NN-amplitudes stem from the isoscalar ci ; 2,3-vertex in the last diagram of 
Fig .2 and they read (modulo polynomials), 
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Note that there are no isoscalar spin-spin and tensor NN-amplitudes Vs,t{q)- This follows from 
the fact that to first order in the 1/M-expansion the Ci^^-vertex remains spin-independent. For 
the sake of completeness we give also the spectral representation of the spin-orbit potential in 
coordinate space (accompanying the standard spin-orbit operator — f (<?i + o-i) • if x V)) which 
reads, 

Vso(r) = -ig f" dfifie-^(l + fir)lmV so ^) . (18) 

With the help of the formula ImL(i/x) = —irk/pi = — (7r/2/i)w/! 2 — 4m 2 the numerical evalu- 
ation of the coordinate space NN-potentials following from the one-loop results eqs. (10-17) is 
straightforward. 

In summary, we have calculated in this work the chiral 27r-exchange NN-potentials at two- 
loop order. We find that these two-loop diagrams lead to a sizeable isoscalar central repulsion. 
Furthermore, we have evaluated here the one-loop 27r-exchange diagram with two second order 
c i,2,3,4- v ertices as well as the first relativistic 1/M-correction to the diagrams with one such ver- 
tex. The analytical results presented here are in a form such that they can be easily implemented 
in a N 3 LO calculation of the two-nucleon system or in an empirical analysis of low-energy elastic 
NN-scattering. 
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